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We construct a non-perturbative method to calculate the effective potential of the \(j>^ theory 
' at finite temperature. We express the derivative of the effective potential with respect to the mass 

, square in terms of the full propagator. We reduce it to the partial differential equation for the 

effective potential using an approximation. We numerically solve it and obtain the effective potential 
non-perturbatively. We find that the phase transition is second order as it should be. We determine 
several critical exponents. 
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, 1 Introduction 

X 

5— ( I It is often expected that broken symmetries are restored at high temperature The temperature- 

^ induced phase transition will be observed in relativistic heavy ion collisions, interior of neutron stars, 

and the early stage of the universe. We may probe new physics through the phase transition at high 
temperature. 

It is, however, very difficult to examine the phase transition. For example, the perturbation theory 
often breaks down at high temperature. As is well-known in finite temperature field theories higher order 
contributions of the loop expansion are enhanced for Bose fields by many interactions in the thermal 
bath |§, H). In the theory physical quantities are expanded in terms of /m^ and XT/m at 
finite temperature. The ordinary loop expansion is improved by resumming the daisy diagram which 

includes all the higher order contributions of O ((ATVm^)") g |, |, |, 1 1, 0, |ll|, ^ |l| The loop 

expansion parameter is XT/m after the resummation. It means that the perturbation theory breaks down 
at T > m/X 1^. Around the critical temperature the ratio m/T is always of 0{X) so a non-perturbative 
analysis is necessary to study the phase transition in Xcj)^ theory|^. 

A variety of methods is used to investigate the phase transition, for example, lattice simulation 
|l||l6[|l^, ^|l§, C.J.T. method H, e-expansion ||, effective three dimensional theory 
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f'l-S'i-H' HI' E3' equation method ||2^, non-purturbative renormalization group method [p9| |3C| , 
|34 |33|] and so on. All the same we still need another method to study the phase transition since they 
are applicable to limited situations. 

In Ref. a new non-perturbative approach was suggested to avoid the infrared divergence which 
appears in the pressure . They differentiated the generating functional with respect to the mass square 
and found the infrared finite expression for the pressure in thermal equilibrium. 

In the present paper we employ the idea of Ref . and develop a new method to calculate the effective 
potential. Differentiating the effective potential with respect to the mass square, we express the derivative 
in terms of the full propagator. We construct the partial differential equation for the effective potential by 
approximating the full propagator. We calculate the effective potential beyond the perturbation theory 
by solving this equation. 

In section 2 we consider the A(/)^ theory at finite temperature and show the exact expression of the 

dV 

derivative of the effective potential — — ^ . We approximate it and obtain the partial differential equation 

for the effective potential. We give the reasonable initial condition to solve this equation. In Sec. 3 we 
solve it and get the effective potential numerically. We obtain the susceptibility, field expectation value, 
and specific heat from it. We determine the several critical exponents by observing their behaviours as 
T varies. The Sec. 4 is devoted to the concluding remarks. 



2 Evolution equation for the effective potential 

As mentioned in Secj^, the loop expansion loses its validity at high temperature. We need a non- 
perturbative method to calculate the effective potential. The effective potential, in general, satisfies the 
following relation, 

dV 

Once we know — — 5- and y(M^), we can calculate the effective potential for arbitrary m^. Following the 

idea, we construct an evolution equation for the effective potential of the Xcj)* theory at finite temperature. 

dV 

In the following we give — — ^ and an appropriate initial condition V(M ). 
We consider the Xcf)^ theory which is defined by the Lagrangian density 

= -I {^y i^^^^' - r"''^' - + + 

where Cct represents the counter term and J is an external source function. If m? is negative, the scalar 
field (j) develops the non- vanishing field expectation value at T = 0. It is expected that the field expectation 
value decreases as T increases and the phase transition takes place at the critical temperature Tc- We 
can explore properties of this phase transition by studying the effective potential at finite temperature. 

Following the standard procedure of dealing with the Matsubara Green function js^, we introduce 
the temperature to the theory. The generating functional at finite temperature is given by 



D[(t>]c^p I dr (fxCs . (3) 




, dV 
In the X(j) theory the derivative of the effective potential — — 5- is expressed by the full propagator of the 

scalar field (See Appendix A), 

dV dVtree BVl dV2 dVct 



dm? dm"^ dm^ dm^ dm^ 



(4) 
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Vtree IS the tree part, 



dVt 



tree ^ ,2 



dm?' 



The non-perturbative effects are contained in Vi and V2, 



dVi 
d'w? 

dV2 



1 



dun? 

Vet is the counter term part 
dV,t 



2^ / 



1 



(27r)3 -p2 _|_ 



|02 +nePo/r _ 1 



dpo 



(27r)3 -p2+p2 



(5) 

(6) 
(7) 



dm? 



{ZmZ^ — 1) 



2 <.+ioo+e 



(27r)3 -p2 + p2 _^ ^2 _^ A^2 + n ePo/'T - 1 
d^p 1 



(27r)3 -p2+p2. 



|02 



n 



(8) 



ay. 



Here 11 = n(p , — pg, 0, m , T) describes the full self-energy. The third term -;- — - on the right hand side 

^ ^ am"' 

of Eq.(^) is divergent. This divergence is removed by the counter term (g) after the usual renormalization 

procedure is adopted at T = 0. The counter term which is determined at T = removes the ultra-violet 

divergence even in finite temperature [0, |l2[ 

We give the initial condition at 0{T^) where the loop expansion is valid, {XT /M A ^ 1). 

We calculate V (M^) by the perturbation theory up to the one loop order. After the renormalization with 

MS scheme at the renormalization scale /2, the one-loop effective potential is 



V{M^) = VtreeiM^) + Vi{M^) + V2{M^) + Vct{M^), 

where Vtree , Vi and V2 + Vet are given by 



VtreeiM^) = ^M^^ 



A 

4!' 



V2iM^)+VetiM^) 



drr^ log 
(M2 H 



1 — exp 



647r2 



log 



(9) 

(10) 
(11) 
(12) 



Note that we need not resum the daisy diagram which has only a negligible contribution of 0{X) for 

In order to investigate the temperature-induced phase transition we consider the theory with non- 
vanishing field expectation value at T = (i.e. m^ takes a negative value, = — ^^). We calculate 
V{—fi'^) with the effective potential (0) by 



y(-M^) = 



For ^ the contribution from 



dm 



' f dVtree ^ 9V ^ d}^ ^ dVet\ 

\ dm'^ 9m2 dm"^ dm^ J 

+ Vtree{M^) + Vi{M^) -f V2{M'') + Vet{M^). 



(13) 



dVi 



is enhanced by the Bose factor. The contribution from Vi can 



be the same order as that from the tree part around the critical temperature. 
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The quantity V2 + will have a negligible contribution: 



fdV2 ^ dV,t\ 



V2(-M') + Kt(V)- 



M2 



drr? + V2{M^) + Vct{M^) 



(14) 



We can show that V2(— /^^) + K:t(— ^^) is really small at the leading order of the loop expansion. At one 
loop level with daisy diagram resummation we find 



647r2 



n 



log 



(15) 



The self-energy satisfies 11 ~ /i^ around the critical temperature for the second-order or the weakly first- 
order phase transition Because we are interested in the effective potential at small (f> region only to 
investigate the phase structure, we neglect ( p^ in the following calculations. 

Furthermore, we ignore the momentum dependence of the self-energy 11 and generate the two point 
function from the effective potential V . We replace as follows in Eq.(Q), 



vn? + -02 + n(0, 0, 0, m\ T)^-^. 



(16) 



We obtain the partial differential equation for the effective potential by integrating over po and angle 
variables in Eq.(^), 



dV 
dm? 



1 

4^ 



drr 



/ 1 



9^ 



(17) 



3 Numerical results 

We calculate the effective potential by solving the partial differential equation (|l^) with the initial con- 
dition Vtree + Vi in Eq. (^ . We solve the equation numerically and show the phase structure of 
theory. 



3.1 analytic continuation 

The integral in (O) is well defined in the region where — =— is real and positive. The effective potential 
^ — ' ocj)^ 

V{<j)) is, however, complex for small <f> below the critical temperature, T < Tc- We have to find the 

analytic continuation in order to calculate the effective potential there. 

To make the analytic continuation, we change the variable of integration r to z through 



j;^ + Z^-Z, (18) 



and rewrite the differential equation (H 



Here Z is the double valued function which is given hy Z — 



9V r,,y^g^. (19) 



1 
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The imaginary part of the effective potential is interpreted as a decay rate of the unstable state [p8| . 

dV 

It is natural that we assume such an imaginary part is negative. The imaginary part of — — ^ should 

be positive in order that the imaginary part of the effective potential may be negative. We have to 

/ 1 d'^V dV 
select the branch of Z = ^ — - — ^ so that imaginary part of — — - will be positive. We calculate the 

effective potential in this branch and the imaginary part of it is always negative as we will see in the next 
subsection. 



3.2 numerical result 

Putting the initial condition Vtree + Vi in Eq.(|^) at = we numerically solve Eq.(p^) and obtain 
the effective potential at rri^ = —fi^. We use the explicit differencing method |3^. In this subsection we 
show the effective potential and calculate critical exponents. 

We illustrate the behaviour of the effective potential at A = 1 in Fig.|l| (a). The field expectation 
value 4>c is the minimum point of the effective potential. It seems to disappear smoothly at the critical 
temperature. We find that the phase transition is second order as it should be. 

For comparison, in Figs.|l| (b), (b') and (c) we show the effective potential calculated by the per- 
turbation theory at one and two loop order with daisy diagram resummation.|^t the one loop order 
an extremely small gap appears at the critical temperature as is clearly seen in Fig.0 (b'). The phase 
transition is first order at the one loop order. 

This situation is modified at the two loop order. We observe no gap and find that the phase transition 
is second order as shown in Figjl] (c). Though Fig.|l| (a) and Fig.^ (c) show the similar behaviour, it will 
be accident. The effective potential calculated up to the two loop order includes the contribution from 
the graphs, Fig.H (a) and Fig.^ (b), with daisy resummation. On the other hand we can take into account 
the contribution from all the other graphs in addition to Fig.^ (a) and Fig.^ (b) within the approximation 
dl^ ) by solving Eq.(p^ automatically. The Fig.|l| (a) accidentally coincides with Fig.|l| (c). 

For T < Tc the effective potential develops a non-vanishing imaginary part at small range. We show 
it in FigJ|. It should be noted that the sign of the imaginary part is always negative. It is consistent 
with the discussion in the previous subsection. 

Evaluating the effective potential with varying the temperature, T, and the coupling constant. A, we 
obtain the critical temperature as a function of A where the field expectation value disappears. We show 
the phase boundary on T-A plain in Fig.^. 

The critical exponents are defined for the second-order phase transition. Around the critical temper- 
ature we expect that the susceptibility x, the expectation value (/)c, and the specific heat C behave as 




Xoc |i|-'',0coc |t|+'%C(x l^r", (20) 

where < = (T — Tc)/T. Analysing the effective potential more precisely we can calculate the critical 
exponents 7, f3, and a. The susceptibility x satisfies the following relation, 

^cxp-i, (21) 

where p is the curvature of the effective potential at 

Since the specific heat C is given by the second derivative of the effective potential around the critical 
temperature, the effective potential V{(l)c) behaves as 

V{^c) oc \tf-". (22) 

We examine the behaviour of p, 0c ,and V{(j)c) around the critical temperature and find the critical 
exponents 7,/? and a. In Figj^ the critical behaviour of p,(l)c and V{<pc) are shown as a function of 
the temperature. We numerically calculate the critical exponents from them. Our numerical results 
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Figure 1: The behaviour of the effective potential V is shown for fixed A(= 1) as the function of the 
temperature. We find no qualitative change for other values A(= 0.5,0.1,0.05). We normalise that 
y(0) = 0. 
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(a) (b) 

Figure 2: Two loop diagrams that contribute to the effective potential. 
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Figure 4: Phase boundary 



Table 1: Critical exponents 





our results 


Landau theory 


experimental results 1 41[ 
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1 


1.24 


a 


- 





0.11 



are presented in the table The critical exponents within our approximation are independent of the 
coupling constant A. We note that our results described in the present subsection remain unchanged even 
when we put the initial mass scale = T^/4 or ~ 4r^. 

4 Conclusion 

We constructed the non-perturbative method to investigate the phase structure of At/)^ theory. The 
derivative of the effective potential with respect to mass square was exactly expressed in terms of the full 
propagator at finite temperature. We found the partial differential equation for the effective potential 
with the replacement (|l^). We gave the initial condition by the 1 loop effective potential in the range 
where the perturbation theory is reliable. We numerically solved the partial differential equation and 
obtained the effective potential. 

Though we made the approximation (p^, we could find that the phase transition of theory is 
second order as it should be. Our method is very interesting because it can show the correct order of the 
phase transition. The approximation ([l^) may be fairy good. 

We determined several critical exponents which roughly agree with those of Landau approximation. 
They are, however, rough values because it is very difficult to solve the nonlinear partial differential equa- 
tion ( p^ ) numerically. We need elaborate a numerical study to obtain more accurate critical exponents. 

^ We use the equations in Rcf|^ to draw them. 
Due to the instabiUty in the exphcit differencing method, we can not see the fine structure of the effective potential 
and can get only the rough values of the critical exponents. We need further numerical study to get more precise values. 
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(c) Minimum of the effective potential V{(j)c) 

Figure 5: Critical behaviour of (/)c, p and V{4>c) 
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The main problem of our non-perturbative method is how to improve the approximation to the 
full propagator. We cannot estimate the error from the approximation (jl^). We need improve the 
approximation to the full propagator in order to know the correction to the current result. 

Our method is very promising since it can probe the region where the traditional perturbation theory 
breaks down. 

Acknowledgements 

The authors would like to thank Akira Niegawa and Jiro Arafune for useful discussions. 



A The derivative of the effective potential in terms of the full 
propagator 



The derivative of the effective potential 
we present details of the calculation of 



dV 

dm? 
dV 

dm? 



can be represented by the full propagator. In this appendix 



given m 



Eq.®. 



We consider the Lagrangian density which is defined by 

2 1 1 



I ( dctioV li—,s2 1 2,2 '^0,4, T i. 



4! 



(23) 



where the suffix denotes the bare quantities. 

We adopt the mass-independent renormalization procedure and represent the effective potential as a 
function of renormalized quantities. The renormalization constants Z and renormalized quantities are 
introduced through transformations 



mo 
Ao 
Jo 



m, 

ZxX, 



(24) 



Using these renormalization constants and renormalized quantities, we separate the Lagrangian density 
( [2^ ) into the tree part L\ and the counter term part Let as ^ 

Le = Ci [0] + £et M + [Ji + Jct)0, (25) 
where Ji + Jet = J ■ The laglangian density Ci and Cct are given by 



dr 



1 

2 
1 
2 

-^(^A^0-i; 



-^(V0)^-imV-^<^^ 



dr 



1 



{Z„,Z^ - l)m202 



(26) 



Here we separate the external source J into Ji and Jet, which satisfy the following equations: 

We expand the field (j){x) around the classical background (f>, 

cj){x) = (f> + ri{x), 



(27) 



(28) 
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and then ij^x) satisfies 



{V)J = 0. 



(29) 



In terms of <j) and rjlx) Eq.(26) is rewritten as 



Cl+Jl<j) = -^771^02 _ ^^4 _^ j^^ 



dri 

a7 



-I v ' ' I'Z \ 2 "^TS ^ 

"ol™ +0'?' ]V -^^V 



3! 



(30) 



and 



A 



-(Z,„Z^ - l)m'^' - ^{ZxZi - 1)0^ + J, 



9?/ 



2^2 



A 



-(Z^Zm - l)m ?7 



A 



where K is defined by 



2 'J""*" ■' 3! 

K = -{ZmZ^ - l)m^ - ^(ZxZl - 1)03 + J^. 



4!' 



From Eq.(|2S|), r/{x) satisfies 

{v)k = 0. 

The generating functional T4^T[t/] for connected Green functions is given by 



Zt[J] = 



(31) 
(32) 
(33) 

(34) 



The effective action Tt{(1)) is defined as the Legendre transformation of Wt[J]. In the spacetime with the 
translational invariance the effective potential V{(f>) is proportional to the effective action. The effective 
potential V{(t>) is 

";y(0) = rT(0) = Wt[j] - S^J, (35) 



T 



where fl — / (fix and the new variable (/> is given by 



5J{y) 



Wt[J] — 0(y) = = const. 



(36) 



Substituting Eqs.(pO[) and (31) into Eq.(|25|), we rewrite the generating functional Zt as a functional of 
renormalizcd quantities, 

Zt = expjp [£i(0) + £,t(0) + (Ji + J,t)0]| 

X j D[r,] exp^ ' dr J d^x {C'M + C'.Av] + Kr,) 



exp { - [Cii4>) + CcS) + (Ji + JcS] } Z'iK). 



(37) 
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Taking into account Eqs.(34) and (p5[), we obtain the effective potential from Eq.(p 

T 



V{$) = -[Ci{$) + Ca 



log Z'{K). 



(38) 



In the mass-independent renormahzation procedure the renormahzation constants Z are independent of 
the mass m. We easily differentiate the effective potential V by the mass square and get 



dVi4>) 
dim? 



+ - l)<^2 



T 1 



DM 

1/T 



(39) 



where we use Eq.(p3|) Q. 

For A(/)^ theory the two-point function (77^(0)) in Eq.(p9) is rewritten as 



l/T 



^ E 

n— — c 
00 

^ E 



dr / (fx5^{x)5T{T){'n{x,T)'q{'Q,Q)) 



(27r)3 



dT J d3a;e*(^P+'^"^)(77(a;,r)7?(O,0)) 



(40) 



(2^)3 iul+p^ + m2 + A 



A ^2 



n' 



where ojn = 27rnT due to the periodic boundary condition for Bose fields and 11 — ll{p^, ^Poi 

the full self-energy. Substituting Eq.(^o|) into Eq.(|39|), we express the derivative of the effective potential 

in terms of the full propagator, 



= 2^ + 2 E y 



dV 
dm? 

+ {ZmZ^ — 1) 



d^p 



I12 V- 

+9 E 



)3 cj2 + p2 _^ ^2 ^ A^2 ^ n 



d3p 1 

(27r)3 ijj'l + p^ + ra^ + ^02 _^ n 



(41) 



Using the residue theorem, we convert the frequency sum T to contour integrals. As long as 11 

n— — oo 

has no singularity along the imaginary po axis, Eq.(^l]) naturally separates into a piece which contains a 
Bose factor and a piece which does not |13, , 



dV 
dm? 



d^p 



1 

47rj 



(2^)3_p2+p2 
1 



- m2 + 1^2 + n ePo/'^ - 1 



+ (^m^0 — 1) 



1 

'27ri 



zoo+e 



(27r)3 -p2 4. p2 ^ ^2 _^ A 02 _^ n 

1 1 

(27r)3 + p2 + ^2 + |02 + n ePo/r - 1 
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1 /•+*°° , f (pp 



47ri 7 (27r)3 -pg + p2 _|. ^2 _,_ A^2 _^ n 
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